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ABSTRACT 


This  subroutine  uses  geometric  altitude  either  in  Metric  or 
English  units  to  calculate  the  corresponding  atmospheric  pressure, 
density,  temperature  in  both  Rankine  and  Kelvin,  speed  of  sound, 
and  acceleration  due  to  gravity. 

This  subroutine  will  compute  in  either  English  or  Metric 
units.  A  control  card  is  used  to  make  this  choice.  The  program  is 
written  in  FORTRAN  11.  This  program  may  be  obtained  from  the 
Scientific  Programming  Section;  Program  No.  SP3<£. 
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LIST  OF  SYMBOLS 


b 

C 

s 

G 

So 

H 

Lk 

M 


Subscript  indicating  basic  or  reference  level 
Speed  of  Sound 

Local  acceleration  due  to  gravity 

Effective  acceleration  of  gravity  at  45°32,40", 
gQ  =  9.  80665  m  sec'^ 

Geopotential  altitude 

Molecular-scale  temperature  gradient 
Molecular  Weight 


m 

rn' 

o 

P 

Q 

R* 


Geometric  meter 

Geopotential  meter 

Subscript  indicating  sea  level  value 

Atmospheric  pressure 
GM0 

Constant  p 

3 

Universal  gas  constant  -  8.  31439  x  10  joules 
(UK)  (K^)-1 


Average  radius  of  earth  t  =  6,  37  1 , 239.  9  meters. 
Real  kinetic  temperature 

Molecular  scale  temperature  where  -  (-^  )T 
Altitude  in  geometric  units 


IV 


Atmospheric  density 

Ratio  of  specific  heats  of  air  Cp/Cy  =  1.4 


v 


PROGRAM  SYMBOLS 


Symbols 

Computer 

Symbols 

Units 

Metric 

English 

G 

GRA  V 

Local  acceleration  due  to 
gravity 

meters/ sec^ 

ft/sec 

H 

11 

Geopotential  altitude 

meters 

ft 

ALT 

ALT 

Geometric  altitude  at  which 
properties  are  desired 

meters 

ft 

T 

M 

ATEMK 

T  emperature 

o 

K 

*T»  O 

R 

ATEMR 

T  emperature 

°R 

P 

APNMM 

Atmospheric  pressure 

Newtons/ 
meters  ^ 

lb/ ft  ^ 

DEN 

Atmospheric  density 

K  /meter^ 
g 

Slugs/  ft 

Cs 

SOS 

Speed  o£  sound 

meters/sec 

ft  /sec 

vi 


INTRODUCTION 


Thib  subroutine  provides  atmospheric  density,  pressure, 
temperature,  acceleration  due  to  gravity,  and  geopotential  altitude 
in  either  Metric  or  English  units.  The  density  is  forced  to  zero  at 
an  altitude  of  400,000  feet.  A  card,  referred  to  as  units,  permits 
this  capability.  If  "units"  ^  0,  the  units  used  are  Metric.  If 
"units"  =  0,  the  units  used  are  English.  Caution:  The  units  of  the 
input  must  be  the  same  as  the  desired  output.  CON50  must  be  equal 
to  0  prior  to  entering  the  subroutine  the  first  time. 

The  program  has  been  checked  out  and  has  less  than  0.  5% 
error  where  compared  to  Reference  1. 

Library  Functions  Used: 


EXPF 


ANALYSIS 


The  equations  listed  below  are  extracted  from  Reference  1 


Geopotential  altitude,  H 


Tm  -  (Tm)k*Lm(H-Ht) 


whtn  Lm  /  0 


V 


p.  Q  -  («-  H*\ 


wUn  Lm  =  0 


e 


**Tm 


Zl 


z 


DISCUSSION 


The  computer  program  is  based  on  Reference  1.  The 
functions  used  are  those  commonly  used  in  main  programs  using 
this  subroutine.  Because  of  this  only  about  1400  digitB  of 
memory  are  required. 


RESULTS 

The  results  of  the  subroutine  have  been  compared  to  the 
data  in  Reference  1.  The  error  is  less  than  0.  5%. 


Start 


Enter  500  from  Page  4 


(Page  6) 
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Enter  206  from  Page  6 


FORTRAN  RUN 


C  PROGRAMMED  BY  W  B  WARREN 

C  TRAJECTORY  ATMOSPHERIC  PROPERTIES 

SUBROUTINE  TAP  *GR A V . AL T , A TE MK , APNMM, DEN , SOS , A TEMRa 
IF  *CON50pU,11,13 

11  DIMENSION  AP*9,4b 
CONSO  K  C0N50  £  1. 

DO  12  I  «  1,9 

12  READ  105,AP*i,ln,AP*I,2B,AP*l,3n,AP*l,4n 
105  FORMAT  S4E15.8P 

READ  100, UNITS 
100  F0RMATXE1S.8O 

13  IF  *UNI TSB500, 1000, 500 
1000  ALT  «  ALT/3. 280839B9 

VL  #  VL/3. 28083989 

500  GRAV  *9.80665»*6371239.9/*6 371239.9  £  AL  Tpp»*2 
H  U  ALT*6371239.9/X6371239.9  £  AL  Tn 
00  202  I  n  2,9 
IF  XH  -  APXl , lau203,202,20? 

202  CONTINUE 
I  *  10 

203  J  *  1-1 

ATEMK  #AP*J,2a£1!M-APi;j,  lon*APXj,3n 
IF  XAPX J , 3bp204»20k,204 

204  APNMM*AP*J,4a»tAPTJ,2n/ATEMKn*.*.03416479/AP*J, 3oa 
GO  TO  206 

APNMM  #  APNMM  •  **. 1 020*2. ?05n  /  «3.281>*2pp 

205  APNMM»AP*J,4a/EXPF?;*. 0341647  4p*XH-APXJ, lOD/APtJ,2ao 
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206  OEN  AS*  .0034838394o*APNMM/ATEMKO 
SOS  A20.046333*ATEMK**.5 
lF*UNITSo50i,502,50l 
502  GRAV  A  GRAV* 3. 28083989 
ALT  *  ALT  •  3.2808389 
HAH*  3.2808389 

DEN  A  OEN  *  **2.20S/*3.28083989**3oo/32. 1740 
SOS  A  SOS  •  3.28283989 
501  ATE  MR  A  SATEMK-273. 16o*S9./5.o  £  491.69 
IF  *UNI TSO600.601 ,600 
601  IF  SAL T-400000. 0504, 504. 505 
600  IF  SALT  -121914. 0604, 6 04 .5 05 
505  OEN  A  0. 

SOS  A  0. 

504  RETURN 
ENO 
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